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Technical  Narrative 


Charles  L.  Dolph  and  Alexander  6.  Ramm 


1.  Scientific  Area,  Central  Objectives,  Main  Approach  to  the  Problem. 

Scalar  and  electromagnetic  wave  scattering  is  of  great  interest  in  physics 
and  technology.  In  the  open  system  theory  in  the  framework  of  the  traditional  ap¬ 
proach  one  should  deal  with  operators  with  continuous  spectrum.  From  the  computa¬ 
tional  point  of  view  it  is  much  more  convenient  to  deal  with  operators  with  discrete 
spectrum.  Furthermore  in  many  practical  problems  the  systems  have  resonant  proper¬ 
ties  so  chat  only  one  or  two  terms  in  an  appropriate  representation  of  solutions  to 
scattering  problems  are  of  importance.  These  ideas  were  extensively  studied  in  the 
physics  and  engineering  literature  [2-5].  Engineering  study  was  initiated  by 
Kacenelenbaum  (1968)  [2]  in  the  USSR  and  by  C.  Baum  (1973)  [3]  in  the  USA.  Mathe¬ 
matical  study  was  initiated  by  Ramm  (1969)  [4];  a  review  paper  which  includes  a 
discussion  of  the  SEM  is  Dolph-Scott  [5].  The  ideas  are  used  in  any  scientific 
area  where  wave  scattering  is  under  study.  The  objective  of  the  theory  is  to  give 
a  computational  approach  to  scattering  problems  which  uses  eigenfunction  expansions 
(of  nonclassical  type)  for  those  operators  with  discrete  spectrum.  These  operators 
are  nonself ad joint  in  most  cases,  because  in  open  systems  there  are  some  losses 
(radiation  losses  etc.)  which  lead  to  dissipation.  Therefore,  from  the  mathematical 
point  of  view,  the  scientific  area  is  the  spectral  theory  of  nonselfadjoint  operators 
and  the  objective  is  the  study  of  eigenfunction  expansions  for  nonselfadjoint  opera¬ 
tors.  These  questions  have  been  insufficiently  studied  in  the  literature  [1].  For 
example,  nothing  is  known  about  eigenfunction  expansions  for  the  basic  operator  of 
laser  theory: 

1  ~ 

Af  -  /  expi(x-y)  f(y)dy 

-1 

From  both  the  mathematical  and  physical  points  of  view  a  very  natural  approach  to 
the  above  problems  consists  of  the  use  of  integral  or,  more  generally,  pseudo-dif¬ 
ferential  equations  on  the  boundary  of  the  scatterer.  This  approach  requires  the 
study  of  spectral  properties  of  Che  corresponding  nonselfadjoint  operators. 


2. 


Air  Force  Relevance  of  the  Research. 


While  such  relevance  is  already  clear  from  the  nature  of  the  problems  men¬ 
tioned  in  Section  1,  it  is  worthwhile  to  point  out  that  the  Singularity  Expansion 
Method  (SEM)  and  the  Eigenmode  Expansion  Method  (EEM)  can  both  be  considered  as 
examples  of  the  general  theory  mentioned  in  Section  1.  Also,  the  problem  of  target 
identification  which  consists  of  the  identification  of  the  obstacle  (target)  from  the 
transient  field  scattered  by  this  obstacle,  can  be  treated  within  the  framework  of 
the  general  theory. 

3.  Main  Results. 

The  following  results  were  obtained  during  the  contractual  period: 

1)  A  criterion  for  the  validity  of  the  SEM  and  EEM  was  established. 

2)  The  basic  relation  between  the  SEM  and  the  EEM  was  given  in  that 
the  poles  of  the  Green's  function  of  the  SEM  were  shown  to  be  in 
one-to-one  correspondence  with  the  zeros  of  the  eigenvalues  of 
the  EEM. 

3)  A  constructive  numerical  process  for  determining  the  poles  of  the 
Green's  function  was  developed. 

4)  The  continuous  dependence  on  the  poles  was  established  for  a  class 
of  boundary  perturbations. 

5)  Sufficient  conditions  for  the  root  vectors  in  the  EEM  to  form  a 
Riesz  basis  [3]  and  a  Riesz  basis  with  brackets  [2]  have  been  found. 

6)  Variational  principles  for  the  spectrum  of  compact  nonselfadjolnt 
operators  have  been  established. 

We  do  not  formulate  these  results  in  detail  because  they  were  reported  to 
AFOSR  earlier  ([6]  -  [10])  and  presented  in  the  series  of  talks  listed  in  the 
Appendix. 

4.  Open  Problems,  Directions  for  Future  Research. 

There  are  many  open  problems  in  this  field.  Many  of  them  are  of  immediate 
practical  interest.  We  mention  some  of  the  open  problems. 

1)  To  what  extent  do  the  complex  poles  of  Green's  functions  determine 
the  shape  of  the  body? 
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2)  To  what  extent  do  purely  Imaginary  poles  of  Green's  functions 
determine  the  shape  of  the  body? 

3)  When  Is  it  more  convenient  to  use  the  single  layer  potential  for 
field  representations  rather  than  the  double  layer  potential? 

4)  What  is  the  short  wave  asymptotic  distribution  of  the  complex  poles? 

5)  What  stationary  principles  exist  for  the  complex  poles? 

6)  What  is  the  relation  between  multiplicities  of  the  complex  poles  and 

of  the  zeros  of  the  eigenvalues  of  the  corresponding  integral  equations? 

7)  To  prove  that  the  root  system  of  the  integral  operators  in  scattering 
theory  forms  a  Rlesz  basis  without  brackets. 

8)  Whether  the  complex  poles  of  Green's  functions  for  convex  lossless 
bodies  are  simple? 

9)  Consideration  of  electromagnetic  scattering  problems. 

5.  Suggestions.  I 

i 

It  would  be  interesting  and  useful  to  organize  some  numerical  experiments  j 
(e.g.,  calculation  of  the  complex  poles)  for  a  practical  evaluation  of  some  of  the 
algorithms  obtained  in  the  course  of  this  research.  To  this  end  a  couple  of  graduate 
students  and  a  postdoctoral  fellow  could  be  incorporated. 
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Appendix  II 


Invited  Talk:  "Nonself adjoint  operators  in  diffraction  and  scattering" 
by  C.  L.  Dolph  and  A.  G.  Ramm;  presented  at  Meeting  of 
American  Mathematical  Society,  Boulder,  Colorado,  March  1980 


NONS EL FADJO I NT  OPERATORS  IN  DIFFRACTION  AND  SCATTERING 


C.  L.  Dolph  and  A.  G.  Ramm 

Department  of  Mathematics 
University  of  Michigan 
Ann  Arbor,  MI  48109 

1.  In  many  problems  of  interest  for  physicists  and  engineers  non- 
selfadjoint  operators  arise  naturally.  For  example  the  EEM 
(eigenmode  expansion  method)  can  be  described  as  follows.  Let  £1 
be  an  exterior  domain  with  a  smooth  closed  compact  boundary  T , 

(A  ♦  k2)u  =0  in  8,  u|r  *  f,  |x|  -  iku)  0  ^ 

2 

as  |x|  -  »  ,  k  >  0  . 

Let  us  look  for  a  solution  of  (1)  in  the  form 

u  *  /fo( t)  exp(ikrxt)  (47rrxt)-1  dt  ,  =  |x  -  t|  .  (2) 

a 

For  the  unknown  a  we  get  equation 

Ao  =  fT  exp(ikrst)  (4  7r  rst)_1  o(t)dt  =  f(x)  .  (3) 

Operator  A  is  nonsel fad  joint  in  L, (F).  Suppose  that  its  root 

2  i 

system  forms  a  basis  of  H  =  1.  (r) .  Then  we  can  look  for  a  solu- 

) 

tion  of  (3)  in  the  form  of  the  scries,  expanding  o(t)  and  f(x) 
according  to  the  root  vectors  of  A(k).  This  is  called  ELM  for 
solution  of  (1).  Thus  we  have  question  1)  Does  the  root  system 
of  A(k)  form  a  basis  of  H  ? 

It  is  easy  to  prove  that  the  Green  function  G(x,y,k)  of  the 
exterior  Dirichlet  (or  Neumann)  problem  can  be  mcromorphic.nlly 


continued  on  the  whole  complex  plane  k  and  its  poles  k  lie 
in  the  lower  half-plane  Im  k  <  0.  If  f  is  a  smooth  function 
with  compact  support,  v  =  /  G(x,y,k) f(y)dy,  /  =  ^  and 
(*)  | vj  <_  C(1  ♦  |kja)_1  ,  a  >  0.5,  C  =  const.,  C  *  C(Im  k),  then 
the  solution  of  the  problem 

Utt  =  Au  in  0  ,  t  >  0,  u|r  =  0,  ujt=0=0,i't  jt,o  =  (4) 

can  be  represented  in  the  form 

u(x,t)=l"=1  exp(-ik^t)Vj(x,t) +o(exp  {- |lm  kn|t}),t  -»•  +«  (5) 

where  v^(x,t)  grow  not  faster  than  tm  as  t  *►  +«  and  m  is 
some  integer.  Expansion  (5)  is  an  example  of  SEM(si ngularity 
expansion  method).  This  leads  to  questions:  2)  What  can  be 
said  about  location  of  the  poles  k^  ?  ;  3)  When  is  (*)  valid?; 

4)  to  what  extent  does  the  set  (k^  }  determine  the  obstacle?;  5)how 
can  the  k^  be  calculated?;  6) whether  the  poles  k^  are  simple? 

2.  The  answer  to  question  1)  is  given  in  [1]  and  is  described 
below.  In  £2] - [4]  some  results  about  bases  with  brackets  are 
given.  Some  answers  to  question  2)  are  given  in  [5] - [ 7] .Answers 
to  question  3)  are  given  in  [8] -[10].  Answer  to  question  4)  is 
unknown.  Answer  to  question  5)  is  given  in  [11], [12].  Answer  to 
6)  is  unknown,  but  some  engineers  (C.E.  Baum  e.g.)  think  that  if 
T  is  convex  then  the  Green  function  of  the  exterior  Dirichlet 
problem  has  simple  poles.  Some  particular  cases  when  this  is 
true  are  discussed  in  [13].  In  [14],  a  survey  of  the  SEM  is  given 
and  [15]  presents  an  engineering  point  of  view.  In  [16]  a  survey 
of  what  is  known  about  questions  l)-6)  is  given  and  in  [17]  some 
relevant  results  can  be  found.  In  [18]  the  relation  between  SEM 
and  the  matchmatical  scattering  theory  is  discussed.  In  [19], [20] 
variational  principles  are  discussed  for  nonscl fad joint  problems. 

• 

3.  In  this  section  we  answer  question  1).  If  A  ■  I,  ♦  T,  where 
L  is  a  sel fad joint  operator  on  a  Hilbert  space  H  ,  the  spectrum 
{A  }  of  L  is  discrete,  A^  *  cn^(l  ♦  0(n  ^))  where  c,p,5  are 


-  7  - 


some  positive  constants,  and  J Tf  J.  <_  |  L  f  | ,  >  0,  a  <  1  for 

all  f  6  D(La),  D(La)CD(T),  then:  1)  the  root  system  of  A 

forms  a  Riesz  basis  of  H  (A  G  R)  if  5  >  1  and  p(l  -  a)  >  2  ; 

2)  if  6  >  0,  p ( 1  -  a)  >  1  it  forms  a  Riesz  basis  of  H  with 

brackets,  (A  6  R^);  3)  if  A*+a(An+j  ®  as 

then  A  6  R.  Let  us  give  the  definition  of  the  Riesz  basis  with 

brackets.  Let  the  system  {f^}  form  an  orthonormal  basis  of  H, 

ml  <  m2  <  '''  mn  <  •••■*■"  is  a  sequence  of  integers,  {F^}  is 

the  sequence  of  the  subspaces,  where  is  the  linear  span  of 

f  ,  f  f  ,  .  Let  { h . }  be  a  minimal  and  complete 

m.  m. +l  m.  ,  -l  1  . 

J  1  J+1 

system  in  H,  H.  is  the  linear  span  of  h»  ,...,h  ,  .  If 

i  m.  m .  ,  ”1  ... 

1  J  +  l 

there  exists  a  map  B  G  L(H) ,  BH^  =  F^,  j  =  1,2, — ,  then  the 
system  {Ik  }  is  called  a  Riesz  basis  of  H  with  brackets  and 
the  numbers  hk  define  the  bracketing.  By  L(H)  we  denote  the 
set  of  linear  bounded  invertible  operators  which  map  H  onto  H. 

4.  Problems. 

—  * 

1.  Find  an  asymptotic  formula  for  k..  as  |Ik|  -*•  «  ,  where 
{kj}  are  the  poles  with  minimal  imaginary  parts.  For  purely 
imaginary  poles  some  information  about  asymptotic  distribution  is 
given  in  [7]. 

2.  To  what  extent  does  the  set  {k^}  determine  the 
obstacle? 

3.  Is  it  true  that  the  complex ‘poles  of  the  Green  function 
of  the  exterior  Dirichlet  problem  are  simple  provided  that  F  is 
convex? 
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Anal.  Appl.  (1980). 
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J.  Math.  Anal.  Appl.  (1980). 

4.  Uniqueness  theorem  of  use  in  general  relativity  (with  P.  Mishnaevsky) , 
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5.  On  the  quasistatic  boundary  value  problem  of  electrodynamics  (with 
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Technical  Narrative 


C.  T.  Tai 


During  the  contract  period  I  have  directed  my  effort  to  the  study 
of  the  singularity  behavior  of  the  eigenfunction  expansions  of  various 
Green's  functions  in  electromagnetic  theory;  they  include  the  Green's 
function  for  Poisson's  equation  and  the  dyadic  Green's  function  for 
Maxwell's  equations.  The  discontinuous  behavior  of  the  eigenfunction 
expansion  of  these  functions  is  discussed  with  the  aid  of  equivalent 
layers  of  surface  charge,  current  sheet,  and  surface  polarization.  The 
principal  result  shows  that  in  an  eigenfunction  of  a  typical  Green's 
function  the  point  singularity  of  a  Green's  function  is  represented  equi¬ 
valently  by  a  layer  of  surface  singularity.  This  characteristic  is 
analogous  to  the  Gibbs'  phenomenon  where  the  representation  of  a  discon¬ 
tinuous  function  by  an  orthogonal  expansion  (Fourier  series)  creates 
the  spikes  which  is  absent  in  the  original  function.  The  only  difference 
is  that  in  our  case  the  surface  discontinuity  has  replaced  the  point 
singularity  in  the  eigenfunction  expansion. 

In  the  case  of  the  dyadic  Green's  function  we  have  examined  in 
detail  the  function  pertaining  to  a  rectangular  waveguide.  Similar  studies 
are  currently  being  done  for  functions  associated  with  other  diffracting 
bodies. 

A  paper,  "Equivalent  Layers  of  Surface  Charge,  Current  Sheet  and 
Polarization  in  the  Eigenfunction  Expansions  of  Green's  Functions  in  Electro¬ 
magnetic  Theory",  based  on  this  research  has  been  written.  It  has  been 
submitted  to  the  IEEE  Transaction  on  Antennas  and  Propagation  for  publication. 

It  is  my  hope  that  this  research  will  be  continuously  supported  by 
the  Air  Force  Office  of  Scientific  Research  so  the  singularity  behavior  of 
various  Green's  function  in  electromagnetic  theory  will  be  better  understood. 


A  communication  has  been  written  to  comment  on  the  recent  paper  by 
A.  D.  Yaghjian  ("Electrical  dyadic  Green's  functions  in  the  source  region", 
Proc.  IEEE,  Vol.  68,  No.  2,  pp.  248-263,  February,  1980).  A  copy  of  the  com¬ 
munication,  which  has  been  accepted  by  IEEE  for  publication,  is  also  attached 
to  this  report.  It  seems  clear  that  much  of  the  confusion  created  by  that 
article  is  due  to  the  lack  of  understanding  of  the  singularity  behavior  of 
dyadic  Green's  functions.  For  this  reason,  it  is  desirable  that  further  work 
should  be  done  to  educate  the  public  as  well  as  to  consolidate  the  concept  and 
the  foundation  of  this  discipline. 
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ELECTRIC  DYADIC  GREEN'S  FUNCTIONS  IN. SOURCE  REGION* 


C.  T.  Tai 


Radiation  Laboratory 
The  University  of  Michigan 
/bin  Arbor,  MI  48109 


The  author  of  a  recent  article  [1]  bearing  the  above  title  has 
criticized  severely  the  derivation  and  interpretation  of  the  well 
established  formula 


E(R)  *  iwvt0  P  §e  •  3  dv*  (1) 

obtained  by  the  delta  function  technique.  Using  the  classical  function 
technique  he  has  derived  a  formula  of  the  form 


E(R) 


iwy  11m 

V° 


/ 


VV« 


3  dv' 


l  •  3 


(2) 


where  denotes  what  he  calls  the  conventional  electric  dyadic  Green’s 
e  - 

function  and  C  Is  the  source  dyadic.  There  are  several  ambiguities  in 
his  paper  that  we  would  like  to  call  attention  to. 

For  convenience  the  function  §e  in  (1)  will  be  designated  as  the 
standard  electric  dyadic  Green's  function  as  distinct  from  §£  in  (2). 
Some  comments  on  the  characterization  of  as  described  by  that 
author  will  be  given  later.  If  In  (2)  is  replaced  by  §e  then 
the  equation  represents  the  so-called  regularized  version  of 
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(1)  as  first  discussed  by  Van  Blade!  [2]  who  only  considers  a  principal 

volume,  V. ,  in  the  shape  of  a  sphere.  In  1972  this  correspondent,  with 

0 

help  from  Or.  Olov  Einarsson,  also  studied  the  source  dyadic  for  bodies 
of  different  shape.  Although  the  work  was  not  published  the  information 
was  passed  verbally  to  various  people  including  the  author  of  [1]. 

What  we  have  done  was  to  demonstrate  the  identity 

lim  /C7  x  v  x  5g  -  k2  5fi]  dv  =  !  (3) 

V +o 

5 

for  several  different  shapes,  V.,  including  spheres  and  cylinders. 

<5 

Equation  (3),  of  course,  is  simply  a  statement  of  the  behavior  of 
in  a  neighborhood  of  the  source  point,  and  this  boundary  condition 
is  implied  when  one  writes  the  differential  equation  for  6g  in  the  form 

v  X  v  X  5e  -  k2  Se  *  I«(R  -  R')  (4) 

using  the  delta  function  together  with  the  idemfactor  as  the 
inhomogeneous  part  of  that  differential  equation.  The  derivation  of 
(1)  based  on  the  delta  function  technique  merely  imposes  this 
characterization  in  an  elegant  way  as  first  shown  by  Levine  and 
Schwinger  [3].  The  formula  itself  without  using  the  dyadic  notation 
can  also  be  derived  by  the  method  of  potentials.  It  is  not  clear  to 
this  correspondent  why  the  author  of  [1]  considers  (1)  to  be  an 
ill-conditioned  solution,  particularly  in  view  of  the  existence  theorem 
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for  this  Integral  discussed  by  FI ki oris  [4].  On  the  other  hand. 


if  one  accepts  (2)  as  the  rigorous  solution  with  5°  replaced  by 
then  It  certainly  can  be  written  in  the  form  of  (1)  since,  by 
definition 


11m  J * '  3^  •  3  dv1  *  -  — —  £  •  3  .  (5) 

v«*°  vs  k* 

It  should  be  mentioned  that  for  numerical  calculations  of  E(R)  involving 
a  free  space  Green's  function  the  regularization  has  to  be  applied  to 
finite  cells.  In  this  case  (1)  can  be  separated  into  two  parts: 


m 


iuu. 


/ 


J  dv1 


V 


/  5e  3JV 


(6) 


The  value  of  the  second  integral  in  (6)  now  depends  not  only  on  the 

shape  of  V.  but  also  on  its  size.  In  a  yet  unpublished  work  [5]  this 

aspect  of  the  problem  is  considered  in  detail  but  the  topic  is  outside 

the  scope  of  the  present  note.  What  we  want  to  point  out  is  that  the 

regularization  with  an  infinitesimal  V.  merely  gives  an  interpretation 

0 

of  (1),  and  in  numerical  work  it  Is  the  finitely  regularized  version 
(6)  that  must  be  used. 

One  of  the  most  confusing  aspects  of  [1]  deals  with  the 
characterization  of  5°.  According  to  [1],  5®  satisfies  the  following 
equations: 


(8) 


I 


7  X  2°  «  g°  . 

6  m 

By  eliminating  2°  from  (7)  and  (8),  one  finds 

v  x  7  x  2°  -  k22°  »  L<s(R  -  R')  (9) 

The  above  equations,  of  course,  is  contradictory  to  (4),  and  the 

solution  to  (9)  is  not  what  the  author  of  [1]  used  in  deriving  (2), 

nor  in  illustrating  Its  application  to  the  waveguide  problem.  It  is 

well  known  that  a  Green's  function,  scalar,  vector,  or  dyadic,  is  an 

entire  entity.  There  is  nothing  like  conventional  or  unconventional. 

It  has  a  unique  and  pre-specified  behavior  in  the  neighborhood  of 

the  singularity,  and  therefore  automatically  satisfies  a  condition  such 

as  the  one  described  by  (3).  It  is,  perhaps,  the  mishandling  of  this 

aspect  of  the  theory  that  leads  the  author  of  [1]  to  create  an 

unorthodox  'equation'  like  (7). 

We  consider  now  the  application  of  (1)  and  (2),  with  2° 

replaced  by  2  in  the  latter  formula,  to  find  E(R)  inside  a  waveguide. 

%» 

The  explicit  expression  for  2g  Is  known  to  have  the  following  form  [6] 

*  S“  -  —  zz<$(R  -  R' )  (10) 

k2 

where  the  series  ?+  applies  to  z  >  z',  the  series  2"  to  z  <  z',  and 
the  source  point  is  located  in  the  plane  z  ■  z'.  The  series  5*  Is 

-  16  - 


v 


discontinuous  at  z  *  z'  and  it  is  also  singular  at  R  ■  R1.  We 
assume  a  current  distribution  which  has  a  z  component  only. 
Substituting  (10)  into  (1)  we  obtain 


E(B)  ■  1u)u0  J  h  ■  3  dv'  ♦  0^  .  (11) 

VJ 

In  evaluating  the  volume  integral  we  use  §+  for  z  >  z'  and  S"  for 
z  <  z'.  The  author  of  [1]  illustrates  the  application  of  (2)  using 
a  principal  volume  in  the  shape  of  a  needle.  He  uses  our  1*  as  his 
5°.  The  result  yields 


E(R)  a  i(uu  Lim 


I 


S±  •  3  dv* 


V°  V,-V, 


a  lorn 

17 

5-  •  3  dv*  -  Lim 

f  S±  .  3  dv1 

0 

J 

L  V, 

V°  " 

(12) 


The  integrals  are  not  carried  out  in  the  referenced  paper,  but  in  a 
private  coranunication  to  this  correspondent  Professor  Robert  E.  Collin 
has  proved  analytically  that  the  second  integral  in  (12)  is  indeed 
equal  to  Jzz/(iuie0)  for  a  typical  distribution  of  Jz  in  a  rectangular 
waveguide.  The  proof  for  an  arbitrary  3  can  be  done  accordingly,  though 
It  is  by  no  means  a  simple  exercise.  In  fact,  this  correspondent 
failed  to  recognize  the  identity  until  it  was  point  out 
to  him  by  Professor  Collin.  From  this  example  it  is  clear  that  there 
Is  no  necessity  to  use  the  regularized  version  of  (1)  to  find  E(R)  in 
a  waveguide.  Unlike  free  space  problems  a  direct  application  of  (1) 


and  (10)  avoids  the  tedious  calculation  of  the  explicit  term  Jzz/(lue0) 
when  the  shape  of  the  principal  volume  Is  not  wisely  chosen.  Of 
course.  If  a  principal  volume  in  the  shape  of  a  pill-box  Is  chosen 
(2)  also  yields  Immediately  (11).  It  should  be  emphasized  again  that 
one  must  use  the  standard  5@  In  (2)  and  not  the  one  described  by 
the  solution  of  (9). 

There  are  many  other  statements  in  [1]  which  are  not  comprehensible 
to  this  correspondent.  For  example,  there  is  absolutely  no  discrepancy 
between  the  work  of  Tai  and  Ro  enfeld  and  that  of  Rahmat-Sami i . 

The  author  appreciates  the  many  valuable  discussions  on  this 
subject  with  Professor  Robert  E.  Collin.  The  kindness  of  Professor  S.  W. 
Lee  for  sending  a  prepublication  copy  of  Reference  [5]  is  gratefully 
acknowledged.  Professor  T.8.A.  Senior  made  valuable  suggestions  in 
regard  to  the  preparation  of  this  communication.  The  comments  from  the 
author's  graduate  students,  particularly  Messrs.  M.  Naor  and  S.  Giles, 
were  very  helpful  in  the  presentation  of  this  work. 
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